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We study ground state properties of the S = 2 quantum antiferromagnetic chain with a bond 
alternation 



3 + 1 



by a Quantum Monte Carlo calculation. We find that the hidden Z2 x Z2 symmetry is broken for 
0.3 < \5\ < 0.5 while it is unbroken in the other regions. This confirms the successive dimerization 
transitions first predicted by Affleck and Haldane. Our result shows that these transitions can 
be understood in terms of the hidden Z2 x Z2 symmetry breaking, as was discussed using the 
Valence-Bond-Solid states. Furthermore, we find that the behavior of the generalized string 
correlation is qualitatively very similar to that in the Valence-Bond-Solid states, including the 
location of zeroes as a function of the angle parameter. 
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Much effpjjt has been devoted to examining Haldane's 



conjectur 



on the difference between quantum anti- 



ferromagnetic chains with integer and half-odd-integer 
spins. These studies also found several interesting phe- 
nomena besides the original conjecture. One of them is 
the successive transitions induced by a bond alternation. 
Consider the Hamiltonian 



(1) 



where 6 represents the strength of the bond-alternation 
(—1 <S<1). The standard spin-S* Heisenberg antiferro- 
magnetic chain corresponds to (5 = 0. When 6 — ±1, this 
Hamiltonian_reduces to that of isolated dimers. Affleck 
and Haldandj-* made an interesting prediction: there are 
25* critical transition points in — 1 < (5 < 1 for spin S. 
(See also ref. ^.) They mapped the model (|^) into an 
0(3) non-linear sigma model with a topological term de- 
fined by the Lagrangian 

C = j-(d^^)id'^^) + ^e^V • X d.^), (2) 

where y> is a three-component vector field with the con- 
straint (f^ ^ 1. The parameters are related to those of 
the original model as follows: 

v = 2SVl~S^, .g=— =L==, a = 2^5(1 + 5), (3) 

where v denotes the spin-wave velocity, which corre- 
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sponds to the "speed of light" in the relativistic field 
theory. Thus the topological angle a depends on the 
bond-alternation parameter S. When S is changed from 
— 1 to 1, a changes from —4ttS to 4ttS. It was argued 
that the above model has a massless (critical) point for 
a = TT (mod 27r), namely 2S critical points are expected 
in model (Q) for -1 < ^ < 1. 

There is another (more intuitive) arguments' based on 
a class of exactly solvable spin chains - Valence-Bond- 
Sohd (VBS) models.i We can introduce VBS states with 
various degrees of dimerization. For S = 2, each spin can 
be constructed by symmetrization of four spin-l/2's and 
we can consider VBS- type states as shown in fig. 1. In 
general, 25* -I- 1 states can be constructed for spin S. If 
we assume that the ground state of the Hamiltonian (|^) 
belongs to the same phase as one of these states, and that 
there are phase transitions between them, there should 
be 25" phase transitions as expected in the field-theory 
argument. 

However, both approaches are based on some assump- 
tions or approximations and thus the prediction should 
be verified by, for example, a numerical calculation. In 
addition, the concept of the order parameter is lack- 
ing in these arguments. How can we distinguish the 
phases which correspond to "different numbers of va- 
lence bonds"? For S = 1, den Nijs and RommelseQ-* 
found that a kind of non-local order parameter charac- 
terizes the Haldane gap phase. (See also ref. ||.) Kennedy 
and TasakiEP revealed that this den Nijs-Rommelse string 
order parameter measures a hidden Z2 x Z2 symmetry 
breaking. One of us extended this idea in ref. ^ to gen- 
eral integer spin, and found that the VBS states can be 
(partly) distinguished by the hidden Z2 x Z2 symmetry 
breaking. Namely, it was found that the symmetry is 
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broken (unbroken) when the number of valence bonds 
is odd (even). Thus the dimerization transitions can be 
understood as hidden Zi x Z2 symmetry breaking, if the 
VBS picture is correct. (See fig. 1.) 

However, the hidden Z2 x Z2 symmetry is not suffi- 
cient for full characterization of all possible VBS states. 
Also it cannot be applied to half-odd-integer spin chains. 
As an attempt to overcome these problems, a general- 
ized string order was proposed.E^r In a dimerized state, 
the generalized string order depends on how the limit is 
taken. Here we choose the definition 



'^string (^) 



lim 



exp (ifii 



2k 

E 

=2j + l 



sns. 



2k+l 



(4) 

This reduces to the Neel order for 6* = 0, and to the den 
Nijs-Rommelse string order for 9 — which is related to 
the hidden Z2 X ^2 symmetry breaking for integer spin. 

It was foundllj) that the generalized string order pa- 
rameter in the (n, m)-VBS state (with n valence bonds 
between sites 2j and 2j + 1 and m between 2j + 1 and 
2j + 2) is given by a polynomial of z = e'^ as 

(5) 



Ostring(^) — /" 



where / is given by 

/'n,m(^) 



2(n + 2)(n+l)^^ 



{2k~~n)z''. (6) 



We note that O"^^.^^^ is not the same for {n,m)- and 
(m, n)-VBS states. The order parameter defined as 
eq. (Q) in the (n, m)-VBS state is equal to an order pa- 
rameter in the (m, rt)-VEiS^ state with a slightly different 
definition. It was provedll2P that the n-th order polyno- 
mial fn,m has n simple roots on the unit circle. Thus the 
generalized string order has n zeroes as a function of 9 
in < < 27r. It was conjectured that this property is 
generic for the partially dimerized phases. We also note 
that 0"ti.ing(^') = <^?tring(2''' " ^) and it is always real. 

For S = 1, the presence of the critical transition 
point |_Lias^ been numerically established by seiieral au- 
thorsMBy^ Recently Yajima and TakahashitJ^ stud- 
ied the S* = 3/2 case and confirmed the presence of the 
transitions. They also measured the generalized string 
order parameter aX 9 — it (this is not related to the hid- 
den symmetry for half-odd-integer spin) and found the 
behavior consistent with the above VBS picture. For 
S ^ 2, the Heiambm-g_^niiiiLt X^r^ 0) was studied by 
several authorsy^BB'BBEl'By) and the behavior of 
the string order parameter was found to be consistent 
with the VBS picture. Tonegawa numerically calculated 
the energy gap of thiL.model (0) with the next-nearest- 
neighbour couplings.E3' 

In this paper, we study the ground state properties of 
the S — 2 chain with the Hamiltonian (|^) more thor- 
oughly by the Quantum Monte Carlo method, with par- 
ticular emphasis on the hidden order and the hidden 
symmetry breaking. 

We peiformed a world-line Quantum Monte Carlo cal- 
culatioiO-' using the Lie-Trotter-Suaiiki product formula 
with checker-board decompositionE^) That is, we made 



an approximation to the partition function Z for tem- 
perature T as Z„ = Tr[(e-'f^^/("^)e-'f^«/("^))"]. Here 
we choose Ha = Ej=odd - = Ej=cvo„ and 
Vj — [1 + d{—iy]Sj ■ Sj+i. The approximate partition 
function Z„ approaches the true partition function Z as 
n ^ 00. We chose a heat-bath algorithm. While we pre- 
pared global fiips along the chain direction, we did not 
use global flips along the Trotter direction and restricted 
the calculation to the ^ = subspace, in order to 
study the ground state. We used the periodic boundary 
condition. For most of the calculation we took T — 0.04 
and n = 96, which turned out to be sufficiently close to 
the extrapolation T and 71 — )■ 00. 

We measured the generalized string correlation (^) at 
9 = mr/20 {n = 0,1,..., 20). We assumed the L de- 
pendence of the correlation function at distance L/2 to 
be 



I L/2-1 \ 

exp (10 J2 ^i )Sl,2 
\ 1=0 I 



A + Bexp(-CL). (7) 



We made a least-square fit of our data for system size 
L{= 40, 60, 80) with the above formula and checked the 
stability of the extrapolation for different choices of the 
system size [L = 40, 50, 60, 70, 80). 

The Heisenberg point: The generalized string order A 
obtained by the fitting (^ and the result (||) for the 

= 2 VBS state are shown in fig. 2. Their qualitative 
behaviors are very similar. In particular, they have ze- 
roes at 9 = 0, TT. The former corresponds to vanishing 
Neel order and the latter means that the hidden Z2 x Z2 
symmetry is unbroken. These results are in agreement 
with the VBS .pictimi.,-ajid, also with previous numerical 
calculationEy^BB-y'BEl'B' at the Heisenberg point. 
The generalized string order has its maximum value at a 
point near 9 — tt/2. The dependence of the generalized 
string-Order on 9 was studied in a short chain by Hat- 
sugaiJlj) The long-range order at = 7r/2 was confirmed 
numerically in refs. ^ and |2^. Here we emphasize that 
the generalized string order is non- vanishing for ^ 0, tt, 
in agreement with the result (||) for the VBS state. We 
observed a similar result at (5 = 0.05. 

The intermediate phase: We show the Neel and the den 
Nijs-Rommelse string correlation functions at (5 = 0.4 in 
fig. 3. The Neel correlation decays faster than that at the 
Heisenberg point. Due to the bond-alternation, the den 
Nijs-Rommelse string correlation has an oscillating part. 
Besides the oscillating part, it seems to have a long-range 
order. This means that the hidden Z2 x Z2 symmetry 
is broken. We show the extrapolated results in fig. 4 
together with the result (|^) for the partially dimerized 
VBS state. The result again agrees (qualitatively) with 
that of the VBS state. At 6* = tt (den Nijs-Rommelse 
string order) there is no longer zero and instead we see 
another zero in < < tt. By symmetry, there are three 
zeroes in < < 27r in agreement with the conjecture 
in ref. |l^. It is somewhat surprising that not only the 
number of zeroes, but also the location of the zeroes in 9 
seems to coincide with the VBS result. For S = 0.3 and 
0.5 we obtained similar results to that for S — 0.4. 

The dimerized phase: Next we consider S = 0.6. As 
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shown in fig. 5, both the Neel correlation and the den 
Nijs-Rommelse string correlation decay to zero. The lat- 
ter means that the hidden Zi x Zi symmetry is unbro- 
ken. We show the extrapolated results in fig. 6, together 
with the result for the completely dimerized state. 
Their behaviors agree well. In particular, the locations 
of zeroes seem to agree again. For b = 0.7 and 0.8, we 
observed similar behaviors to that for b — 0.6 (but more 
similar to the result for the completely dimerized state, 
as expected). 

To summarize, we measured the generalized string 
correlation in the S' = 2 spin chain with a bond- 
alternation (|l|) by the Quantum Monte Carlo method. 

We observed that there are three phases in < (5 < 1, 
confirming the prediction by Affleck and Haldane. The 
hidden Z2 x Z2 symmetry is broken only in the inter- 
mediate phase. This is consistent with the VBS picture 
of the transitions. Hence the transitions can be under- 
stood as hidden Z2 x Z2 symmetry breaking. From this 
argument, it would be natural to expect that the tran- 
sitions belong to the same c = 1 universality class as in 
the S = \ case. This argument is again consistent with 
the field-theory prediction by Affleck and Haldane. 

Furthermore, we found that the behaviors of the gen- 
eralized string order in those phases are quite similar to 
those in the VBS-type state. In particular, the locations 
of the zeroes seem to agree. It might be possible to de- 
termine the underlying mechanism for this phenomenon. 
(For example, the zeroes may be related to some un- 
known hidden symmetries.) 

Our result shows that the region < 5 < 0.05 be- 
longs to the S = 2 Haldane phase, 0.3 < 6 < 0.5 to the 
intermediate phase, and 0.6 < (5 < 1 to the dimerized 
phase. Thus the transition points should be located in 
the regions 0.05 < 5 < 0.3 and 0.5 < i5 < 0.6. However, 
in these regions, the correlation length seems to grow 
and hence we could not make reliable extrapolation and 
could not estimate the transition points more precisely. 

From eq. (H) we can obtain the field-theory prediction 
for the transition points: 5 = 0.25 and S — 0.75. The 
latter is clearly inconsistent with our result 0.5 < 6 < 0.6. 
(The value-wasjipt consistent with the numerical result 
for S = IBBEJ and for S = 3/2^) In general, we 
cannot expect that the field theoretical argument gives 
the quantitatively precise location of the critical point. 
However, the qualitative prediction of the occurrence of 
the phase transition and the description of the critical 
behavior are reliable. 
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Fig. 1. The VBS picture of the 5 = 2 successive dimerization 
transitions. Each dot represents a spin 1/2. The bonds represent 
singlet pairs ("valence bonds"), and the broken circle indicates 
the symmetrization of four spin-l/2's to form a spin 2. The 
number of valence bonds between the neighboring sites changes 
as the dimerization proceeds. The hidden Z2 x Z2 symmetry is 
broken only in the intermediate state. 

Fig. 2. The extrapolated generalized string order at (5 = as a 
function of 8. The exact result in the VBS state is also shown. 
Their behaviors are quite similar, and both have 6 = 0, tt as 
zeroes. The relatively large error is presumably due to the large 
correlation length compared to the cases for <5 = 0.4 and 0.6 

Fig. 3. (a) The Neel and (b) the den Nijs-Rommelse string cor- 
relations at (5 = 0.4 for system sizes L = 20, 40, 60 and 80. The 
den Nijs-Rommelse string correlation is long-ranged, while the 
Neel correlation decays exponentially. 

Fig. 4. The extrapolated generalized string order at <5 = 0.4 as a 
function of 0. The exact result for the partially dimerized VBS 
state is also shown. Their behaviors are quite similar, including 
the location of the zero. 

Fig. 5. (a) The Neel and (b) the den Nijs-Rommelse string cor- 
relations at (5 = 0.6 for system sizes L = 20, 40, 60 and 80. Both 
correlation functions decay to zero. 

Fig. 6. The extrapolated generalized string order at 5 = 0.6 as a 
function of 8. The exact result in the completely dimerized state 
is also shown. Their behaviors are quite similar, including the 
location of the zero. 
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